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Abstract. We introduce a new class of nonlocal kinetic equations and nonlocal Fokker- 
Planck equations associated with an effective generalized thermodynamical formalism. These 
equations have a rich physical and mathematical structure that can describe phase transitions 
and blow-up phenomena. On general grounds, our formalism can have applications in different 
domains of physics, astrophysics, hydrodynamics and biology. We find an aesthetic connexion 
between topics (stars, vortices, bacteries,...) which were previously disconnected. The common 
point between these systems is the (attractive) long-range nature of the interactions. 

1. Introduction. The statistical mechanics of systems with long-range interactions 
is currently a topic of active research in physics . Systems with long-range interactions 
are numerous in nature: self-gravitating systems, two-dimensional vortices, non-neutral 
plasmas, metallic clusters, dipoles, fracture etc. These systems exhibit similar features 
such as negative specific heats, inequivalence of statistical ensembles, phase transitions, 
self-organization and persistence of metastable states. Among all the previous examples, 
self-gravitating systems and 2D vortices play a special role because they both interact 
via an unshielded Newtonian potential (in dimensions D = Z ov D = 2) and possess 
a rather similar mathematical structure |2] ■ Coincidentally, the chemotactic aggregation 
of bacteries in biology has some connexions with the collapse and organization of self- 
gravitating systems and 2D vortices E] . 

For systems with long-range interactions, the mean-field approximation turns out to 
be exact in a proper thermodynamic limit [1]. Therefore, their equilibrium description 
amounts to solving a variational problem, namely the maximization of the Boltzmann 
entropy 5*3 [/] at fixed mass M = M[f] and energy E — E[f]. In the case of quantum 
particles, the Boltzmann entropy is replaced by the Fermi-Dirac entropy Spoif]- This 
maximization problem determines the most probable distribution of particles at statisti- 
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cal equilibrium, assuming that all the accessible microstates (with given E and A4) are 
equiprobable. However, if, for some reason, the microstates are not equiprobable, other 
forms of entropy can emerge. This appears in the physics of complex media displaying a 
fractal (or multifractal) structure and exhibiting anomalous diffusion. In general, these 
effects are due to the existence of "hidden constraints" that change the form of the 
transitions probabilities that we would naively expect. Motivated by this problem, we 
shall develop an effective generalized thermodynamical formalism (in /i space) associated 
with a larger class of entropy functionals. Specifically, we consider the maximization of a 
generalized entropy S[f] = — / C{f)d^rd^v, where C is a convex function, at fixed mass 
M and energy E. Since the energy is fixed, we call this a microcanonical description. We 
also consider the minimization of a generalized free energy F = E — TS at fixed mass 
M and temperature T. Since the temperature is fixed, we call this a canonical descrip- 
tion. We discuss the inequivalence of these two descriptions when the caloric curve (3{E) 
presents turning points or bifurcations. This occurs in particular when the potential of 
interaction is long-ranged. 

We also introduce a new class of relaxation equations associated with this generalized 
thermodynamical formalism [3 0. We first consider a generalized class of Fokker- 
Planck equations extending the ordinary Kramers and Smoluchowski equations JSj . These 
equations have a canonical structure as they decrease a generalized free energy. When 
the potential of interaction is long-ranged, these equations are non-local in space and can 
exhibit a rich variety of behaviors including phase transitions and blow-up phenomena. 
In the limit of short-ranged interactions, they reduce to Cahn-Hilliard equations [7]. We 
also consider a generalized class of kinetic equations extending the ordinary Boltzmann 
and Landau equations . These kinetic equations have a microcanonical structure as 
they conserve energy and increase a generalized entropy. 

Our formalism can have applications in different domains of physics, astrophysics, 
fluid mechanics, biology, economy etc. with various interpretations that are not necessar- 
ily connected to thermodynamics. It is therefore important to develop a general formalism 
without reference to a specific context. Then, a justification has to be given in each case. 
For example, the maximization of the functional S[f] = — J C{f )dPrdP'v, where C is 
convex, at fixed mass and energy determines a nonlinearly dynamically stable stationary 
solution of the Vlasov-Poisson system for coUisionless stellar systems (and 2D vortices). 
In this context, iS'[/] is a H-function 8 , not an entropy. Then, our generalized relaxation 
equations B| can be used as numerical algorithms to construct stable stationary solu- 
tions of the Vlasov-Poisson system. On the other hand, equations similar to generalized 
Fokker-Planck equations appear in biology in relation with the chemotactic aggregation 
of bacterial populations. In any case, it is useful to develop a thermodynamical analogy 
and use a vocabulary borrowed from thermodynamics. Thus, we can directly transpose 
the methods developed in thermodynamics to a different context. 

2. Maximum entropy principle Let us consider a system of N particles in interac- 
tion and denote by /(r, v, t) their distribution function defined such that fd^rd^v gives 
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the mass of particles with position r and velocity v at time t. The spatial density is 

(1) p{r,t) = J f{r,v,t) d^v, 
and the total mass 

(2) " / 

Let F(r,t) = — V$ be the force (by unit of mass) experienced by a particle. We 
assume that the potential <i>(r, t) is related to the density p{r, t) by a relation of the form 



(3) 'i>(r,t)= / p{v',t)u{r~v')d 



where u{r — r') is an arbitrary binary potential depending only on the absolute distance 
|r — r'l between the particles. The energy can be expressed as 

(4) E = J ^fv^d^rd^v +^Jp^d^r = K + W, 

where K is the kinetic energy and W the potential energy. The following results remain 
valid if <i> = ^ext (r) is a fixed external potential, in which case the potential energy reads 

We introduce a generalized entropy of the form 

(5) S ^ - J C{f) d'^rd^v, 

where C(/) is a convex function, i.e. C"{f) > 0. We are interested by the distribution 
function /(r,v) which maximizes the generalized entropy (0 at fixed mass and energy, 
i.e. 

(6) Max S[f] I E[f] = E, M[f] = M. 

Since the energy is fixed, we shall associate this maximization problem to a microcanonical 
description. Introducing Lagrange multipliers and writing the variational principle in the 
form 

(7) 6S - I35E - a6M = 0, 

we find that the critical points of entropy at fixed mass and energy are given by 

(8) C'(/) = -/3e-a, 

where e = + is the energy of a particle by unit of mass. The Lagrange multipliers 
(3 = 1/T and a can be interpreted as a generalized inverse temperature and a generalized 
chemical potential, respectively. Equation (jSJ can be written equivalently as 

(9) / = F(/3e + a), 
where F{x) — (C")~^(— a;). From the identity 

(10) /'(e) = -/3/C"(/), 

resulting from Eq. (jSJ, we find that /(e) is a monotonically decreasing function of energy 
if /3 > 0. The case of negative temperatures /3 < can also be of interest depending 
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on the form of the function C . Exphcating the relation between the potential and the 
density, the equilibrium distribution is determined by the integro-differential equation 

(11) = j /(r',v')^(r-r')d^r'd^v'|-a. 

The conservation of angular momentum L = J fr x v(Pr(fiv (in D — 3) can be easily 
included in the variational principle (O by introducing an appropriate Lagrange multiplier 
pn, where $1 is the angular velocity. Equation (jSJ remains valid provided that e is replaced 
by the Jacobi energy ej — e — H ■ (rxv) — i + f where w = v — Jlxris the relative 
velocity and <&e// = $ — |(Oxr)^ is the effective potential accounting for inertial forces. 
We now introduce the generalized free energy 

(12) j[f]^s[f]-mf], 

associated with the functional and © . We are interested by the distribution function 
/(r,v) which maximizes the generalized free energy (jl2|l at fixed mass and temperature, 
i.e. 

(13) Max J[f] ^ S[f] ~ mn I M[f]^M. 

Since the temperature is given, we shall associate this maximization problem to a canon- 
ical description. Introducing Lagrange multipliers and writing the variational principle in 
the form 

(14) SJ - a6M = 0, 

we find that the critical points of free energy at fixed mass and temperature are given by 

(15) C'{f)^-Pe-a, 
as in the microcanonical description. 

3. Stability conditions The critical points of the variational problems © and (|13(l 
are the same. The equilibrium state is then obtained by solving the integro-differential 
equation l(TT|l and relating the temperature (3 to the energy E. We can thus plot the 
generalized caloric curve (3{E) parameterizing the series of equilibria. We need now to 
determine the stability of these solutions by investigating the sign of the second order 
variations of S or J. In the microcanonical situation, we must select maxima of S[f] at 
fixed mass and energy. The condition that / is a maximum of S at fixed mass and energy 
is equivalent to the condition that 5^ J = 5^S — (36'^ E is negative for all perturbations 
that conserve mass and energy to first order. This condition can be written 

<5V = ~ 1 C"(/)^d^rd^v - i/?! <5p<5<i>d^r < 0, 

(16) V (5/ I 6E = SM = 0. 

The condition of stability in the canonical situation requires that / is a maximum of J[f] 
at fixed mass and temperature. This is equivalent to the condition that S^J is negative 
for all perturbations that conserve mass. This can be written 

5V = - y C"(/)^d^rd^v - l/?! SpS^d^'r < 0, 
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(17) V 5f I SM = 0. 

Using Eq. ()10f) . the functional arising in these stability criteria can be expressed as 

(18) d\j - i/?! / iffid^rd^v - / 5p<5<i>d^r 



r U /'(e) 

We emphasize the importance of the first order constraints in the stability analysis. We 
note that canonical stability imphes microcanonical stabiUty but that the converse is 
wrong in general. Indeed, if inequality H17|) is satisfied for all perturbations that conserve 
mass, it is a fortiori satisfied for perturbations that conserve mass and energy. Since the 
converse is wrong, this implies that some solutions can be stable in the microcanonical 
ensemble while they are unstable in the canonical one. The microcanonical and canonical 
descriptions are inequivalent when the caloric curve f3{E) presents turning points. This 
situation is well-known in the case of systems with long-range interactions such as self- 
gravitating systems |0 . The stability of the solutions can be decided by using the turning 
point criterion of Katz jlOj which extends the theory of Poincare on linear series of 
equilibria. It is found that a change of stability in the series of equilibria occurs in the 
microcanonical ensemble when the energy is extremum and in the canonical ensemble 
when the temperature is extremum. Stability is lost or gained depending on whether 
the series of equilibria P{E) turns clockwise or anti-clockwise at that critical point. An 
illustration of these results is proposed in Ref. 0, in the case of self-gravitating fermions. 
A change of stability along a series of equilibria can also occur at a branching point llOj , 
where the solutions bifurcate. A general classification of phase transitions for systems 
with long-range interactions has been proposed by Bouchet & Barre |11| . 

4. The free energy functional The maximization problem Ijl^^ll in the canonical 

ensemble can be simplified. First of all, we write the free energy in the usual form 

(19) F[f]=E[f]-TS[f]. 
To solve the minimization problem 

(20) Min F[f]^E[f]^TS[f] \ M[f] = M, 

we shall proceed in two steps. We first minimize F[f] at fixed density p(r). Introducing 
a Lagrange multiplier A(r), we find that the global minimum /, (r,v) of this variational 
problem is determined by 

(21) ^?'(/*) = -/3 y + A(r) 
The distribution function can be written 

(22) /. 

where F{x) = {C')^^{—x). We define the density and the pressure by 



/3( y + A(r^ 



(23) p = / /d^v, p=^ I fv^d^w. 
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Substituting Eq. H22(l in the foregoing expressions, we find that 

(24) p = ^5(/3A), p = -^KPX), 
with 

r-\-oc 

(25) g{x)=2^SD f F{x + t)t^dt, 

Jq 

1 r+oo 

(26) h{x)^—2^SD F{x + t)t^dt, 

D Jo 

where Sd is the surface of a unit sphere in D dimensions. Ehminating A between the 
foregoing expressions, we find that the fluid is barotropic, in the sense that p — p{p) 
where the equation of state is entirely specified by C{f). We can now express the free 
energy (|19|l as a functional of p by using F[p] = F[f^]. The energy Q is simply given by 

(27) E^^ j pd^v+]^j p^d^v. 
On the other hand, the entropy (0 can be written 

(28) S ^ -^jjj^ J d^r C[F{t + /3A)] dt. 
Integrating by parts and using C'[F{xy\ = —x, we find that 

(29) S = J d^'r F'{t + /3A)(t + pX)t''/^dt. 

Integrating by parts one more time and using Eqs. (|24|l . (|25|1 and H26() . we finally obtain 

(30) S = ^^^P J Pd^r + P J Xpd^v. 
Collecting all the previous results, the free energy (|19|l becomes 

(31) F[p] = ^ I p(^X+^y^r+^l p^d^r. 

Finally, using the relation h'{x) — —g{x) obtained from Eqs. H25() and H26() by a simple 
integration by parts, it is easy to check that Eq. H24I) implies 

(32) x+'--- r^-^dp'. 



Hence, the free energy can be written more explicitly as 

(33) = / ^ ^dp'd^r + \ I p^d^v. 
We are led therefore to study the minimization problem 

(34) Min F[p] I M[p] = M, 

for the free energy functional H33|) . Writing the first order variations in the form 

(35) 5F - aSM = 0, 

we find that the critical points of free energy satisfy the condition of hydrostatic balance 

(36) Vp = -/9V$, 
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and that A = $ + a//?. We can easily check that the condition of hydrostatic equihbrium 
is directly implied by the relation / = /(e) derived in Sec. 2. This immediately results 
from the identities 

(38) ^p^hj ^"^"^""^ = hj z'w^'f^"'^''^ = -p^^- 

Explicating the relation between the potential and the density, the condition of hydro- 
static equilibrium can be written in the form of an integro-differential equation 

(39) ^^Vp - -V / p(r')w(r - r')d^r'. 

P J 

This equation determines the equilibrium solutions (critical points of F). Their stability 
(minima of F) is determined by the condition 

(40) V dp I SM = 0. 

Using the condition of hydrostatic equilibrium, the above functional can be rewritten as 



(41) ^^^-27^^ ^^W^^^"^^ 

5. Examples of entropy functionals Among all functionals of the form (O, some 
have been discussed in detail in the literature. The most famous functional is the Boltz- 
mann entropy 

(42) SB[f] = - I flnfd^'rd^v. 

It leads to the isothermal distribution 

(43) / = Ae-P^ 

The corresponding distribution is physical space is the Boltzmann distribution 

/r,^N D/2 

(44) p = A'e-^*, A. 

The distribution function H43() leads to the classical equation of state 

(45) p = pT, 
and to the free energy 

(46) F[p]^T J p\np d^r+^ J p<^> d^r. 

Closely related to the Boltzmann entropy is the Fermi-Dirac entropy 

(47) SMf] = - f{-lr,l+(l- I] Infl - 1) L^rd^v, 

which leads to the Fermi-Dirac distribution function 

(48) / = 
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The Fermi-Dirac distribution function 1)48(1 satisfies the constraint / < ?7o which is related 
to Pauli's exclusion principle in quantum mechanics. The isothermal distribution function 
((4^{|l is recovered in the non-degenerate limit / ^ 770- The distribution in physical space 
associated with the Fermi-Dirac statistics is 



(49) P= ""^ff. I^-A^e^% 
where /„ is the Fermi integral 

(50) Ut)=l ^d.. 

The distribution function ((48|l leads to the quantum equation of state given in parametric 
form as 

(51) 



P 



pD/2 



-In i(A), p= n IniA). 



Recently, there was a considerable interest in physics for functionals of the form 
(52) S,[f]^--^J{r-f)d^rd^^, 

where g is a real number. Such functionals introduced by Tsallis jl2| are called g-entropies. 
They lead to polytropic distribution functions 



(53) 



/ 



(g-1)/? ■ 

/i e 



The index n of the polytrope is related to the parameter q by the relation n ^ D/2 + 
l/{q — 1) . Isothermal distribution functions are recovered in the limit g — s- 1 (i.e. 
n +00). The distribution function (|53|l leads to the polytropic equation of state 



(54) 
with 
(55) 



p = Kp'' 



1 

n 



K = 



-l/n 



where A = [{q — l)/3/g]5-i and B{n,m) is the Beta function. The corresponding free 
energy can be written 

(56) F[p] = -^^j[p-'-p)d^v+\j d^r, 
and the density in physical space is 

(57) p 



7-1 
A- i $ 



We note that a polytropic distribution with index q in phase space yields a polytropic 
distribution with index 7=1-1- 2(g — l)/[2 -I- D{q — 1)] in physical space. In this sense, 
Tsallis distributions are stable laws. By comparing Eqs. (|53l) and H57|) or Eqs. (|19|l . (|52|l 
and (|56|) we note that K plays the same role in physical space as the temperature T — 1/ (j 
in phase space. It is sometimes called a "polytropic temperature" . 
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We have started a systematic study of the variational problems (jSJ , H13|l and (|34|l by 
considering first the gravitational interaction in D dimensions. The case of isothermal 
self-gravitating systems associated with the Boltzmann entropy has been considered in 
|14j . The case of self-gravitating fermions associated with the Fermi-Dirac entropy has 
been considered in [HIEI- The case of polytropic self-gravitating systems associated with 
the Tsallis entropy has been considered in fP?) . 

6. Physical applications The maximization problems discussed previously can have 
various physical applications (not necessarily related to thermodynamics) that we briefly 
mention. 

(i) Statistical mechanics: The variational principles JHl, and 1)34(1 determine 
the statistical equilibrium states of systems with long-range interactions, such as self- 
gravitating systems. In that case, S[f] is the Boltzmann entropy (|42|l for classical particles 
(e.g., stars in globular clusters) ^Hl or the Fermi-Dirac entropy (|47|l for fermions (e.g., 
massive neutrinos in dark matter models) [J. The microcanonical situation applies 
to isolated Hamiltonian systems such as stellar systems ^7] . The canonical situation (|13|) 
or l|34|) applies to systems in contact with a heat bath imposing its temperature, like for 
the interstellar medium This is also the proper statistical description of a gas of 
self-gravitating Brownian particles EH] . This discussion remains valid for other types 
of long-range interactions. 

(ii) Vlasov equation: The variational principle © determines nonlinearly dynamically 
stable stationary solutions of the Vlasov-Poisson system 

(59) A$ = AttG J fd^v, 

which describes collisionless stellar systems. These robust organized states can emerge 
as a result of a violent relaxation '22[ I23|. In that context, S[f] is called a H- function 
[Sj. Boltzmann and Tsallis functionals are particular H- functions (not true entropies in 
that context) associated with isothermal and polytropic stellar systems Note that 

the criterion of nonlinear dynamical stability l(H|l- H16() is more refined than the stability 
criterion ((13(l - ((17|l . These results on nonlinear dynamical stability remain valid when the 
force F is related to the distribution function / by a relation of the general form Q . 

(Hi) Euler- Jeans equation: The variational principle (|34|) determines the nonlinear 
dynamical stability of stationary solutions of the Euler-Jeans-Poisson system 

(60) 1^ + V • (pu) ^ 0, 

(61) -^-^(u- V)u = --Vp-V$, A$=47rG/9, 

describing barotropic stars. In that context, yV[p] — F[p\ + J p^cPr represents the energy 
of the star. The first term in Eq. (|33|1 is the internal energy of a barotropic gas and the 
second term is the gravitational energy (a kinetic term has also to be introduced in 
the energy functional W). Since the stability condition (|34|l implies ©, we note that a 
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stellar system is stable with respect to the Vlasov equation whenever the corresponding 
barotropic star with the same equation of state is stable with respect to the Euler- Jeans 
equations (but the reciprocal is wrong in general) |2H . In astrophysics, this is known as 
the Antonov's first law |24j . 

(iv) Generalized thermodynamics: The variational principles l|13|l and (|34|l deter- 
mine the generalized thermodynamical stability of complex systems exhibiting anomalous 
diffusion. Such systems have a complicated phase space structure (fractal, multifractal,...) 
due to the action of microscopic constraints (hidden constraints) that are often difficult to 
formalize. These systems can be described by effective equations resembling generalized 
Boltzmann and Fokker-Planck equations f3! IB] . The variational principles , H13|l and 
(|34|l also determine the linear dynamical stability of stationary solutions of these equa- 
tions. Accordingly, there is a close connexion between thermodynamical and dynamical 
stability 0- 



7. Generalized Kramers equation We shall introduce formally a relaxation equa- 
tion associated with the minimization problem (|2()|l . We write this equation in the form 

df 



(62) 



dt ' dr dv dv 



where S/Sf denotes a functional derivative. By construction, this equation conserves mass 
since the right hand side is the divergence of a current —J/- Using the expression of the 
free energy (|19|l with Eqs. 10} and lO, we obtain the generalized Fokker-Planck equation 



(63) 



D 



/C"(/)|^+/3/v 



where L is the advection operator. Morphologically, Eq. (|63|l can be viewed as a gener- 
alized non-local Kramers equation j25| . It describes the dynamics (in phase space) of a 
system of Langevin particles in interaction governed by a generalized class of stochastic 
processes 



V,f/(ri,...,rjv) + W2i?/, 



R,(t), 



where — ^^v^ is a friction force and [/(ri, rjv) — J2i<j'^i^i ^ ^j) ^ potential of 
interaction. The last term is a generalized stochastic force. The usual white noise term 
R^(i), satisfying {Ri{t)) = and (i?a,i(t)i?6,i(t')> = S^J6abSit - t'), where a,6 = 1, ...,D 
refer to the coordinates of space and i,j = 1, ...,N to the particles, is multiplied by a 
function that can depend on the local (macroscopic) distribution function fi = fijct, "Vijt). 
When C(/) = / In/, which is related to the Boltzmann entropy (|42|1 . Eq. H64(l reduces 
to the usual Langevin equations 

(65) ^--^v, -V,;i7(ri,...,rjv) + \/iDR.(i). 



dt " dt 

This stochastic process describes a system of Brownian particles in interaction. More 
generally, the term in front of Ri(i) in Eq. (|64|l can be interpreted as a multiplicative 
noise since it depends on r, v. Note that it depends on r, v through the distribution 
function /(r,v,t). Therefore, there is a back-reaction from the macrodynamics, leading 
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to a situation of anomalous diffusion. In that context, the generalized Kramers equation 
can be derived from the N-body Fokker-Planck equation by using a Kramers-Moyal 
expansion and a mean- field approximation 0113. The first term in Eq. H63|l is a generalized 
diffusion (depending on the distribution function) and the second term is a friction. The 
friction coefficient ^ = D(3 satisfies a generalized Einstein relation. Note that D can 
depend on r, v and t without altering the general properties of the equation. We can use 
this indetermination to write the generalized Kramers equation in the alternative form 

d^^ C"{f) 

which will have the same general properties as Eq. H63() . This equation involves an ordi- 
nary diffusion and a nonlinear friction. Equation (|66|) can be deduced from Eq. H63() by 
the substitution D' = DfC"{f). One of these two forms will be prefered depending on 
the situation contemplated. 

The generalized Kramers equation H63|l can also be obtained from a variational prin- 
ciple by maximizing the rate of free energy dissipation F at fixed mass and temperature 
[3]. Therefore, Eq. (|63|l satisfies a canonical H-thcorem F < Q, provided that _D > 0. 
Indeed, an explicit calculation yields 



(66) 



(67) 



/C"(/)§^+/3/v 
av 



d^vd^w < 0. 



DT 

~r 

This shows that F is the Lyapunov functional of the generalized Kramers equation. Now, 
at equilibriimr F = 0, so that according to Eq. H67(l . 

(68) ^^+/3v = 0. 

av 

Integrating with respect to v, we get 

(69) C"(/) = -/3^+A(r). 

The cancellation of the advective term Lf = in Eq. (|63|) combined with Eq. (|69|l implies 
that / = /(e) and VA = — /3V$. Therefore, A{r) = — /3$(r) — a and we recover Eq. lfTK|l . 
This shows that a stationary solution of Eq. H63() extremizes the free energy F at fixed 
mass and temperature. In addition, only minima of F at fixed M and T are linearly 
stable with respect to the generalized Kramers equation (|63|l . Indeed, considering the 
linear stability of a stationary solution of Eq. H63() , we can derive the general relation 

(70) 2Xd^F = 6^F < 0, 

connecting the growth rate A of the perturbation Sf ~ e-^* to the second order variations 
of the free energy F and the second order variations of the rate of free energy production 
5^F < 0. Since the product X5^F is negative, we conclude that a stationary solution 
of the generalized Kramers equation H63|l is linearly stable (A < 0) if, and only if, it is a 
minimum of free energy F at fixed mass and temperature. This aesthetic formula shows 
the equivalence between dynamical and thermodynamical stability for the generalized 
Kramers equation. Therefore, it only selects minima of F, not maxima or saddle points. 
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8. Generalized Smoluchowski equation We shall now introduce formally a relax- 
ation equation associated with the minimization problem H34(l . We write this equation in 
the form 



(71) 



dp 

dt 



= V 



1 



f5F_ 



Using the expression of the free energy we obtain the generalized Fokker-Planck 

equation 



(72) 



dp 
dt 



Morphologically, Eq. H72|l can be viewed as a generalized non-local Smoluchowski equation 
155 . It can be obtained from a variational principle by maximizing the rate of free energy 
dissipation F at fixed mass and temperature 3 . Therefore, Eq. H72|l satisfies a canonical 
H-theorem F < 0. Indeed, an explicit calculation yields 



(73) 



pV$)2d^r < 



This shows that F is the Lyapunov functional of the generalized Smoluchowski equation. 
Now, at equilibrium _F = 0, and we recover the condition of hydrostatic balance H3t)|) . It 
is also possible to show 3 that only minima of free energy F are linearly stable with 
respect to the generalized Smoluchowski equation (|72|l . 

Explicating the relation between the potential and the density, the generalized Smolu- 
chowski equation H72|l can be written 



(74) 



dp 
dt 



V 



p'{p)\/p + p\/ / u{r-r')p{r',t)d 



The Lyapunov functional associated with this equation is the free energy 
p{p') 



(75) 



F[p\ = j pj^ ?^dp'd^v+\ j p{v,t)u{v 



D„jD„f 



r')p{r',t)d^rd 



The stationary states are determined by the integro-differential equation 

P'ip), 



(76) 



Vp = -V / p(r')u(r-r')d^r'. 



Equation (|74|l generalizes the usual Smoluchowski equation in two respects. First, the 
coefficient of diffusion ^p'{p) is an arbitrary function of the density p associated with 
a generalized entropy functional (first term in Eq. (|75|l '). Second, this equation is non- 
local due to the presence of an arbitrary binary potential of interaction u(r — r') in the 
energy functional (second term in Eq. H75() '). This equation possesses therefore a very rich 
mathematical and physical structure. 

We have started a systematic study of the generalized Smoluchowski equation lf71|l by 
considering first the gravitational interaction in D dimensions. The case of an isothermal 
equation of state corresponding to the Boltzmann free energy has been considered in 
\W\ 1141 1261 n7\ . It describes a system of self-gravitating Brownian particles (see a review 
in [201) ■ The case of a quantum equation of state associated with the Fermi-Dirac free 
energy has been considered in |28[ 0] . It describes a system of self-gravitating Brownian 
fcrmions. Finally the case of a polytropic equation of state associated with Tsallis free 
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energy has been considered in jl3j. It describes a system of self-gravitating Langevin 
particles experiencing anomalous diffusion. 



9. The method of moments The Smoluchowski equation can be derived from 
the Kramers equation in a high friction limit ^ +00, or equivalently for large times 
t 3> £,~^- This can be shown easily by using a method of moments. Integrating Eq. H63(l 
over velocity, we get the continuity equation 

(77) |^+V.(pu)=0, 

where u = (1/p) / fv(Pv is the local velocity. Multiplying Eq. H63() by v and integrating 
over velocity, we get the momentum equation 



d d 9 9$ 

(78) ^{pui) + - — {pu^Uj) + - — Pij + p- 



D 



fC"{f)^+Pfv, 



dt dxj dxj dxi 

where Pij — J fwiWjCpv is the pressure tensor and w = v — u the relative velocity. 
Introducing the notation (f){f) = xC" {x)dx, the first term in the right hand side can 
be rewritten d(j){f)/d'v and, since it is a gradient of a function, it vanishes by integration. 
We are left therefore with 



(79) 



d_ 

dt 



9$ 

^dxi 



-S.pui 



In the high friction limit ^ = Df3 — > -|-cxd, assuming [3 of order unity, the term in 
bracket in Eq. H63|l must vanish so that the distribution function satisfies 

,,2 



(80) 



C'{f) = -(3 



+ A(r,t) 



where A(r,t) is related to the spatial density p{r,t) through the relation 
(81) p = [ fd^w. 



P - 



We find therefore that u = 0{£, ^) and Pij = p6ij + 0(^ ^) where p{r,t) is the local 
pressure 

D 

determined by Eq. (|80|l . As in Sec. 4, the fluid is barotropic, i.e. p{r,t) — p[p{r,t)]. To 
first order in the momentum equation H79|l implies that 

1 



(82) 



(83) 



pu 



:(yp + p\7<P) + 0{C 



Inserting the relation (|83|l in the continuity equation H77|l , we get the generalized Smolu- 
chowski equation H72|l . The generalized Smoluchowski equation, as well as the first order 
correction to the distribution function /(r, v, i), can also be obtained from a formal 
Chapman-Enskog expansion \7\ . 



10. GeneraUzed kinetic equations The Kramers equation and the Smoluchowski 
equation have a canonical structure in which the temperature is fixed. We shall now 
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introduce a generalized kinetic equation possessing a microcanonical structure in which 
the energy is fixed jSI E] ■ This is the generahzed Landau equation 



(84) 

with 
(85) 



dt 



^v^' 



u 



where u = vi — v is the relative velocity and ^ is a constant. We can also consider the 
alternative form 



(86) 



The generalized Landau equation can be derived from a generalized Boltzmann equation 
in a weak deflexion approximation. In turn, the generalized Kramers equation H63|l can be 
derived from the generalized Landau equation in a thermal bath approximation jH]. The 
generalized Landau equation satisfies a H-theorem (S* > 0) for the generalized entropy 
The entropy © plays therefore the role of a Lyapunov functional. At equilibrium, 
S = 0, and we obtain the distribution (jS)). In addition, it can be shown jH| that only 
maxima of S at fixed M and E are linearly stable with respect to the generalized Landau 
equation. 

Generalized kinetic equations such as H84(l appear when the transition probabilities 
are different from the form that we would naively expect due to the action of microscopic 
constraints (hidden constraints). One particular case is when the particles are fermions. 
In that case, the "hidden constraints" correspond to the Pauli exclusion principle that 
prevents two particles with equal spin to occupy the same phase space cell. This exclusion 
has been explained by quantum mechanics. More generally, the "hidden constraints" may 
not have necessarily a fundamental origin. 



11. Particular examples It can be of interest to discuss some special cases explicitly. 
For the Boltzmann entropy H42|l . C"(/) = 1// and Eq. 1)63(1 reduces to the ordinary 
Kramers equation 



(87) 



a/ 

dt 



d_ 
dv 



D 



dl 
dv 



-/3/v 



The corresponding equation in physical space, obtained in the high friction limit, is the 
ordinary Smoluchowski equation 



at 



i(rvp + /3V$) 



Finally, Eq. ((84|l leads to the ordinary Landau equation 



(89) 



dl 
dt 



Lf- 



d 



df ^dh 



For the Fermi-Dirac entropy H47|) . C"{f) = i/fiva ^ /)• Iii order to avoid the diver- 
gence of the term fC"{f) as / ^ 770, it is appropriate to consider the alternative form 
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of the generalized Kramers equation. This yields 



(90) 



5/ 
dt 



d_ 
dv 



D' 



a/ 

9v 



+ fifim - /)v 



The corresponding equation in physical space is given by Eq. (|72|l with the equation of 
state (|51f) . In fact, when the Kramers equation is written in the form Ht)()|) . the coefficient 
in factor of the diffusion current in Eq. (|72|l is not 1/^ but is more complex. We refer to 
[7] for a detailed discussion of this subtle point. Finally, Eq. H8t)|l leads to the fermionic 
Landau equation 



(91) 



dl 
dt 



dvi^ 



- fiVo - /) 



dv'^ 



For the Tsallis entropy (|5^ . C"{f) = qf ^ and Eq. has the form of a nonhnear 
Kramers equation 



(92) 



dl 
dt 



Lf = 



d_ 
dv 



D 



dfl 
dv 



The corresponding equation in physical space is the nonlinear Smoluchowski equation 



(93) 



dt 



i(ii:vp'' + pv$) 



Finally, Eq. H84|l leads to the (/-Landau equation 



(94) 



dl 
dt 



Lf = 



d 
dvi^ 



dp dfi 



dv'- 



dvl 



12. Maximum entropy principle in physical space We now consider a system 
of N particles in interaction with total mass 



(95) 

and energy 
(96) 



M ^ j p{r,t) d^r, 



where the potential <i>(r, t) is related to the density p(r, t) by a relation of the form 
(97) $(r, t)^ J pir', t)u{r - r')d"r'. 



We introduce a generalized entropy in physical space 



(98) 



S 



C[p) d^r, 



where C{p) is a convex function, i.e. C"{p) > 0. We are interested by the distribution 
p(r) which maximizes the generalized entropy H98(l at fixed mass and energy, i.e. 

(99) Max S[p] I E[p] = E, M[p] = M. 

Since the energy is fixed, we shall associate this maximization problem to a microcanonical 
description. Introducing Lagrange multipliers and writing the variational principle in the 
form 

(100) SS - f3SE - aSM = 0, 
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we find that the critical points of entropy at fixed mass and energy are given by 

(101) C"(p) = -/3$ - a. 

The Lagrange multipliers /3 = l/T and a can be interpreted as a generalized inverse 
temperature and a generalized chemical potential, respectively. Equation (|101(l can be 
written equivalently as 

(102) p = F(/3$ + a), 
where F{x) — {C')^^{~x). From the identity 

(103) /,'($) = ~(ilC"{p\ 

resulting from Eq. H1Q1|I . /o(<i>) is a monotonically decreasing function of $ if /? > and 
a monotonically increasing function of $ if /3 < 0. Explicating the relation between the 
potential and the density, we find that the equilibrium distribution is determined by the 
integro-differential equation 

(104) C'{p) = -/3 / p(r')w(r - v')(fv' - a. 



We now introduce the generalized free energy 

(105) J[p]^S[p]-l3E[pl 

associated with the functionals H9()|l and H98|) . We are interested by the density p(r) which 
maximizes the generalized free energy (|105|) at fixed mass and temperature, i.e. 

(106) Max J[p\ = S[p] - l3E[p] \ M[p] = M. 

Since the temperature is given, we shall associate this maximization problem to a canon- 
ical description. Introducing Lagrange multipliers and writing the variational principle in 
the form 

(107) 6 J - a6M = 0, 

we find that the critical points of free energy at fixed mass and temperature are given by 

(108) C"(p) = -/3$ - a, 
as in the microcanonical description. 

13. Stability conditions In the microcanonical situation, we must select maxima 
of S[p] at fixed mass and energy. The condition that p is a maximum of S at fixed 
mass and energy is equivalent to the condition that S^J = S'^S — is negative for all 

perturbations that conserve mass and energy to first order. This condition can be written 

<5V = - 1 C"(p)^d^r j <5pJ<&d^r < 0, 

(109) y 5p\ 5E = 5M = Q. 

The condition of stability in the canonical situation requires that p is a maximum of J[p\ 
at fixed mass and temperature. This is equivalent to the condition that 5^ J is negative 
for all perturbations that conserve mass. This can be written 

- f C"(p)^d^r - i/3 / Spd^fd^'r < 0, 
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(110) V 5p I 5M = 0. 

Using Eq. H1(J3|) . the functional arising in these stabihty criteria can be expressed as 



(111) 



The discussion on the inequivalence of the microcanonical and canonical descriptions is 
the same as the one given in Sec. 3. 

14. Examples of entropy functionals The Boltzmann entropy in physical space 
(112) Sb[p\ = - I plnpct'^Y, 



leads to the isothermal distribution 
(113) p = Ae~l"^. 

The Fermi-Dirac entropy in physical space 

leads to the Fermi-Dirac distribution 



This distribution satisfies the constraint p < ag which puts an upper bound on the local 
density of particles. It usually arises when one considers finite size effects or when one 
introduces a lattice model in physical space preventing two particles to occupy the same 
site. The isothermal distribution pi3(l is recovered in the dilute limit p ^ ctq. 
The Tsallis entropy 

1 

9~1 

where g is a real number, leads to the polytropic distribution 

1 

The isothermal distribution (|113|l is recovered for g ^ 1, i.e. n 



(116) /(p?_p)d^r, 



(117) 



15. Physical applications The maximization problems discussed previously can 
have various physical applications (not necessarily related to thermodynamics) that we 
briefly mention. 

(i) Statistical mechanics of point vortices: The variational problems H99|l and I|1U6I) 
arise in the statistical mechanics of 2D point vortices provided that we make the corre- 
spondance p ^ uj between the density and the vorticity and the correspondance $<-!■'(/' 
between the potential and the stream function. In that context, the maximization problem 
(|99|l associated with the Boltzmann entropy H112|l determines the statistical equilibrium 
state (most probable distribution) of a cloud of point vortices [221 Q] ■ 
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(ii) Two-dimensional Euler equation: The maximization problem 199(1 determines non- 
linearly dynamically stable stationary solutions of the 2D Euler-Poisson system 

(118) -^ + u-Vcj = 0, 

(119) AV-^-w, 

where u = — z x Vi/' is the velocity field. The 2D Eulcr equation describes the evolution 
of the density distribution of point vortices in the "coUisionless" regime (Vlasov regime) 
before correlations (collisions) have developed (501 . It also describes the inviscid evolution 
of continuous vorticity flows in two-dimensional hydrodynamics. In that context, S[lij\ is 
called a H-function. Boltzmann and Tsallis functionals are particular H-functions asso- 
ciated with isothermal and polytropic vortices Note that the criterion of nonlinear 
dynamical stability H1U9|) is more refined that the stability criterion (|110|) which is itself 
more refined than the Arnol'd theorems |^ 1221 0] • 

(Hi) Generalized thermodynamics: The variational principles ((99|l and ((l()6|l determine 
the generalized thermodynamical stability of complex systems exhibiting anomalous dif- 
fusion in physical space. Anomalous diffusion can have various origins: it can can be 
due to the existence of traps, of a lattice preventing particles to reach occupied sites, or 
any other constraints that are often difficult to formalize. This is what we call "hidden 
constraints" ^ . These systems can be described by effective equations resembling gener- 
alized Fokker-Planck equations [H]. The variational principles ((99|l and p06|l determine 
the linear dynamical stability of stationary solutions of these equations 



16. Generalized drift-difFusion equation We shall introduce formally a relaxation 
equation associated with the maximization problem HlUt)|) . We write this equation in the 
form 

dp 



(120) 



where 5/5p denotes the functional derivative. By construction, this equation conserves 
mass since the right hand side is the divergence of a current — Jp. Using the expression 
of the free energy ((l()5|l with Eqs. ((96|l and ((98|l . we obtain the generalized drift-diffusion 
equation 

dp 



(121) 



dt 



D 



pC"{p)Vp + (3pV<^ 



Morphologically, this equation can be viewed as a generalized non-local Smoluchowski 
equation. It describes the dynamics (in physical space) of a system of Langevin particles 
in interaction governed by a generalized class of stochastic processes ^^: 



(122) 




IDp, 



C{pi 



R,(t), 



where pi = p{ri,t) and Ri(t) is a white noise satisfying {'Ri{t)) — and {Ra.i{t)Rb,j{t')) — 
5ij5ab5{t — t'), where a,b — 1, ...,D refer to the coordinates of space and i,j = 1, ...,iV 
to the particles. The particles interact via the potential C/(ri, ...,rjv) — ^Ylii<j ~ '^j) 
and C{p) is an arbitrary convex function. When Ci^p) = phip, which is related to the 
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Boltzmann entropy H112(l . Eq. (|122|l reduces to the usual Langevin equations 
(123) ^ = -eV.C/(ri,...,rA.) + ViDR,(i). 

This stochastic process describes a system of Brownian particles in interaction. The term 
in front of Tli{t) in Eq. (|122|l can be interpreted as a multiplicative noise since it depends 
on the position r. Note that it depends on r through the density p{r, t). Therefore, there is 
a back-reaction from the macrodynamics. In that context, the generalized Smoluchowski 
equation H121|l can be obtained from the A^-body Fokker-Planck equation by using a 
Kramers-Moyal expansion and a meanfield approximation 0] . The first term in Eq. 
(I121|l is a generalized diffusion (depending on the density) and the second term is a drift. 
The drift coefficient ^ — D(3 satisfies a generalized Einstein relation. Note that D can 
depend on r and t without altering the general properties of the equation. We can use 
this indetermination to write the generalized drift-diffusion equation in the alternative 
form 



(124) 



dt 



D' 



-V4> 



C"{p) 

which will have the same general properties as Eq. (|121|1 . This equation involves an 
ordinary diffusion and a nonlinear drift. Equation (|124|) can be deduced from Eq. H121() 
by the substitution D' = DpC"{p). One of these two forms will be prefered depending 
on the situation contemplated. 

The generalized drift-diffusion equation \121\ can also be obtained from a variational 
principle by maximizing the rate of free energy production J at fixed mass and tem- 
perature Therefore, Eq. H121() satisfies a canonical H-theorem J > 0. An explicit 
calculation yields 



(125) 



J= - 



pC"{p)Vp + (3pV^ 



d^r > 0. 



If P > 0, this inequality is true provided that _D > 0. If p can take both signs (which 
depends on the initial conditions), it is more convenient to consider Eq. (|124ll . In that 
case, the inequality J > is true provided that D' > 0. This shows that J is the Lyapunov 
functional associated with the generalized drift-diffusion equation H121() . 
Now, at equilibrium J = 0, so that according to Eq. (|125|l . 

(126) VC"(p) + /3V$ = 0. 
Integrating, we get 

(127) C'{p) = -/3$ - a, 

and we recover Eq. 1)108(1 . Therefore, a stationary solution of Eq. ((121|l extremizes the free 
energy J at fixed mass and temperature. In addition, it is shown in j2j that only maxima 
of J at fixed M and T are linearly stable with respect to the generalized Fokker-Planck 
equation H121(l . 

Explicating the relation between the potential and the density, the generalized drift- 
diffusion equation ((124|l can be written 



(128) 



dp 
dt 



D 



Vp 



/3 



C"(p) 



M(r- r')p(r',i)d 
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with the free energy 

(129) J = - J Cip) d'^r-^P J p<^> d^r. 

This equation generahzes many drift-diffusion equations introduced in the hterature (we 
refer to [6'6\ for a connexion with mathematical works and for a detailed list of references) . 
It would be of interest to investigate its properties by remaining as general as possible, 
i.e. without specifying the function C{p) and the binary potential of interaction u(r — r') 
explicitly. For example, which conditions must satisfy C{p) and u{r — r') to generate 
blow-up solutions? Can we regroup the functionals S[p] is "classes of equivalence" with 
the heuristic idea that functionals of the same "class" will yield "similar" results? Which 
functionals S[p] lead to confined solutions with a compact support, such as polytropic dis- 
tributions associated with the Tsallis functional 23|? These are interesting mathematical 
problems which could be tackled in relation with Eq. H128() . 

17. Relation to Cahn-Hilliard equations. If we now consider the case of short- 
range interactions, it is possible to expand the potential 

(130) ^{r,t) = J u{r')p{r + r')d°r' 
in Taylor series for r' ^ 0. Introducing the notations 

(131) a = Ju{\x\)d"x and b = ^ J u{\x\)x^d^x , 
we obtain to second order 

(132) <^{r,t) = ap{r,t) + ^Ap{r,t). 
In that limit, the free energy takes the form 

(133) J[p]^lpb J + 

where we have set V{p) = —{2/bf3)C{p) — {a/b)p'^. This is the usual expression of the 
Landau free energy. In general, (3b is negative so that we have to minimize the functional 
integral. For systems with short-range interactions, the conservative equation H128(l be- 
comes 

(134) |^ = v||pV(Ap-y'(p)) 

This is the Cahn-Hilliard equation which has been extensively studied in the theory of 
phase ordering kinetics. Its stationary solutions describe "domain walls" . We can view 
therefore Eq. H128() as a generalization of the Cahn-Hilliard equation to the case of systems 
with long-range interactions. Therefore, its general study is of great mathematical and 
physical interest. 



18. Physical applications The generalized Fokker-Planck equations discussed pre- 
viously and in Pl El can have various physical applications (not necessarily related to 
thermodynamics) that we briefly mention. 
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(i) Numerical algorithms: The relaxation equations presented in Sees. 7, 8 and 16, 
and in Refs. JH!, can be used as numerical algorithms to solve the variational problems 
(0), p3(l . H34(l ■ H99|l and (|l()6|l . This is of great practical interest because it is in general 
difhcult to solve the integrodifferential equations l|ll|l . (|39|l and p()4|l directly and be 
sure that the solution is a maximum of S or J . These relaxation equations can thus be 
used as numerical algorithms to construct thermodynamically stable equilibrium states 
of Hamiltonian systems with long-range interactions (see Sec. 6 (i) and 15 (i)) as well 
as nonlinearly dynamically stable stationary solutions of the Vlasov, Euler- Jeans and 2D 
Euler equations (see Sees. 6 (ii), (iii) and 15 (ii)). 

(ii) Violent relaxation: The relaxation equations introduced in [341 1231 1^1^ provide a 
small-scale parameterization of the Vlasov-Poisson and 2D Euler-Poisson systems. They 
describe the convergence of the flow, on a coarse-grained scale, towards a metaequilibrium 
state corresponding to a galaxy in astrophysics or a large-scale vortex (e.g., Jupiter's 
great red spot) in 2D hydrodynamics. The theory of violent relaxation is discussed in, 



(iii) Statistical mechanics: In the isothermal case, the nonlocal Kramers and Smolu- 
chowski equations describe the evolution of the distribution function and density of a gas 
of Brownian particles in interaction [201 El- They can be seen, therefore, as the canonical 
counterpart of the nonlocal Vlasov, Boltzmann and Landau equations describing the evo- 
lution of the distribution function of a Hamiltonian system of particles in interaction, for 
which a microcanonical description holds 5 . Since statistical ensembles are not equiva- 
lent for systems with long-range interactions, it is of conceptual interest to compare the 
microcanonical and canonical descriptions to see the analogies and the differences 1191. 

(iii) Biological colonies: Non-local drift-diffusion equations also occur in biology, in 
connexion with the chcmotactic aggregation of bacterial populations. A model of chemo- 
tactic aggregation has been proposed by Keller & Segel . In some approximation [JTj , 
their equations reduce to the Smoluchowski-Poisson system as for self-gravitating Brow- 
nian particles |19[ I20| . This analogy is developed in jj. In a more general context, the 
diffusion coefficient or the chemotactic drift can depend on the density. This can take 
into account anomalous diffusion or finite size effects preventing unphysical blow-up of 
the bacterial concentration. A simple regularization of the Smoluchowski-Poisson system 
is provided by the equation 



which is associated with the Fermi-Dirac entropy 1)1141) in physical space. This equation 
respects the constraint p < a aX any time and therefore prevents blow-up. By rescaling the 
diffusion coefficient D{r, t) appropriately, this equation can be put in the form l|72() with 
an effective pressure p — Tln(l — p/cr). A more general model of chemotaxis is provided 
by the drift-diffusion equation (|128|l which can take into account anomalous diffusion and 
finite size effects of various forms. As we have seen, this equation is associated with an 
effective thermodynamical formalism. An even more general model is represented by the 



e.g., isanniEniEi. 



(135) 
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non-Markovian equation 
(136) I^V.ji. 



/3 







This equation can take into account delay effects that are relevant in the problem of 
chemotaxis (indeed, the equation satisfied by <& is of the form edt^ = A$ — Xp — fj,^ 
|36| ) . We note that Eq. (|136l) does not admit a Lyapunov functional (or a free energy) so 
that the effective thermodynamical formalism developed previously breaks up. It would 
be of interest to see how these ideas can be generalized to that context. 

(iii) Generalized thermodynamics: The relaxation equations presented in Sees. 7, 8 
and 16, and in Refs. can be considered as effective kinetic equations which may be 

useful to model complex systems. As indicated previously, they can arise when the sys- 
tem is subject to microscopic (hidden) constraints that change the form of the transition 
probabilities. The consequence is that the accessible microstates arc not equiprobable, 
resulting in new forms of entropic functionals. An explicit and fundamental example 
is played by the Pauli exclusion principle in quantum mechanics which leads to the 
Fermi-Dirac entropy instead of the Boltzmann entropy. We have also exhibited a gener- 
alized stochastic process (|64ll - (|122ll which leads to generalized Fokker-Planck equations 
of various forms. These stochastic equations form just a particular example of processes 
generating anomalous diffusion and complex phase space structure. More general and 
more realistic microscopic processes could also be considered and studied. However, since 
the generalized kinetic equations JHSJ, l|72(l and (|121|l can be obtained from arguments 
of a very wide scope, such as the Maximum Entropy (Free Energy) Production (Dis- 
sipation) Principle 3^ for example, we believe that they have a relatively fundamental 
and universal character and that they will be obtained from a large class of microscopic 
processes. 



19. Conclusion In this paper, we have developed a generalized thermodynamical for- 
malism and corresponding kinetic theories. We believe that generalized thermodynamics 
is relevant for the physics of complex systems when "hidden constraints" are in action 
[H]. These microscopic constraints imply that the microstates with given mass and energy 
are not equiprobable, contrary to what is usually postulated in thermodynamics. We can 
either use the Boltzmann entropy and try to take into account these additional micro- 
scopic constraints, or consider only the usual constraints (mass and energy) and change 
the form of entropy [^. This introduces some indetermination that is encapsulated in 
the q parameter of Tsallis or more generally in a function C{f). We emphasize that this 
indetermination is intrinsic to the problem and not a flaw of our approach. It occurs 
because we do not have a complete knowledge of the system's dynamics. However, we 
have suggested |3] that generalized entropies can be regrouped in "classes of equivalence" 
and that, for a given system, a class is more appropriate than another. This extends the 
attempt of Tsallis and co-workers to try to find the correct value of q corresponding to a 
given (complex) system. We argue that, more than a scalar g, a full function C{f) must 
be considered in general. 

On the other hand, the maximization problems {Max S \ E,M fixed} and {Min F = 
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E — TS I M fixed} can arise in different situations that are not necessarily related to 
thermodynamics. For example, the maximization of a _ff function at fixed mass and 
energy determines nonlinearly dynamically stable stationary solutions of the Vlasov- 
Poisson system in astrophysics. In that case, this maximization problem is related to 
dynamics, not thermodynamics. However, in order to investigate this dynamical problem, 
it can be useful to develop a thermodynamical analogy and use a similar vocabulary. We 
have given other examples (in fluid mechanics, biology,...) where this thermodynamical 
analogy could be developed. 
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